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ABSTRACT 

We use the most extensive integral field spectroscopic map of a local galaxy, NGC 628, 
combined with gas and stellar mass surface density maps, to study the distribution 
of metals in this galaxy out to 3 effective radii (Ro). At each galactocentric distance, 
we compute the metal budget and thus constrain the mass of metals lost. We find 
that in the disc about 50% of the metals have been lost throughout the lifetime of 
the galaxy. The fraction of metals lost is higher in the bulge (~70%) and decreases 
towards the outer disc (~ 3 Ro). In contrast to studies based on the gas kinematics, 
which are only sensitive to ongoing outflow events, our metal budget analysis enables 
us to infer the average outflow rate during the galaxy lifetime. By using simple physi¬ 
cally motivated models of chemical evolution we can fit the observed metal budget at 
most radii with an average outflow loading factor of order unity, thus clearly demon¬ 
strating the importance of outflows in the evolution of disc galaxies of this mass range 
(log(M*/M 0 ) ~ 10). The observed gas phase metallicity is higher than expected from 
the metal budget and suggests late-time accretion of enriched gas, likely raining onto 
the disc from the metal-enriched halo. 

Key words: galaxies: abundances - galaxies: evolution - galaxies: fundamental pa¬ 
rameters - galaxies: individual (NGC 628) 


1 INTRODUCTION 


Current models of galaxy evolution predict that star for¬ 
mation and consequent metal enrichment in galaxies are 
regulated by a combination of gas inflows and outflows. In 
this framework, the enrichment of the intergalactic medium 
(IGM) is mediated by metal-rich outflows produced by 
galactic winds, especially during the early, more active 


phases of galaxy evolution (Oppenheimer & Dave 2006 
Scannapieco et al.|2006[[Shen, Wadsley fc Stinson|2010||Op-| 


penheimer et al.|2012 I. Evidence for ongoing metal-enriched 


outflows has been obtained both through absorption line 
studies of quasars whose lines of sight cross the vicinity of 
nearby galaxies (Tumlinson et al. 2011[ 2013 Werk et al. 


|2013[ [Bregman et al.||2013[ ) and through direct imaging of 
outflowing ionised or neutral/molecular gas at low and high 


redshift (Heckman et al. 

2000| Martin||2005| Veilleux, Cecil 

& Bland-Hawthorn|2005| 

Martinl20061 fTremonti , Moustakas 

& Diamond-Stanic 2007 

IFeruglio et al.pOlOl Sturm et al. 

20111 Gano-Dfaz et al. 20121 Maiolino et al. 2012 Gicone 


et al.||2015l. However, these observations are only probing 


ongoing events. A global understanding of the effect of out¬ 
flows on the cosmic star formation history and IGM enrich¬ 
ment is only possible if we can gain insight into gas flows 


over the entire evolutionary history of galaxies. The con¬ 
tent and distribution of metals in galaxies can provide such 
information. 


Metals are produced by stellar nucleosynthesis and act 
as a fossil record of a galaxy’s star formation history and gas 
flows. Simple analytical models relating metal abundances 
to star formation rate (SFR) and gas flow rates have been 


presented by several authors ( 

Erb 2008 

Lilly et al. 

2013 

Dayal, Ferrara & Dunlop||2013 

|Peng & k 

aiolino||2014| |As-| 

casibar et al. 

2014 

1 with the aim of explaining the scaling 


relations observed in statistical samples of galaxies. More 
detailed chemical evolutionary models have been success¬ 
fully exploited to reproduce the observed chemical abun¬ 
dances, gas content and SFR of the solar neighbourhood, 
the Milky Way or nearby disc galaxies (Golavitti et al.|2009 


Marcon-Uchida, Matteucci fc Gosta||2010| [Spitoni fc Mat-| 
teucci|2011 1. 


The common aim of chemical evolutionary models is 
to use chemical abundances observed today to infer the gas 
flow history and star formation history of the system. For 
example, in order to reproduce the chemical abundances of 
the Milky Way, a two-stage inflow model is generally in¬ 
voked. The bulge is assumed to have formed at early times 
during a first inflow phase, while the disc is assumed to 
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2 F. Belfiore et al. 


form inside-out during a second inflow phase ( 

Matteucci & 

Francois||1989| [Boissier & Prantzos||1999| [Naa 

3 & Ostriker 

2006 

Williams et al.|2009 

1. Hydrodynamical simulations of 


disc galaxies generally agree with the inside-out disc forma¬ 


tion model (Pilkington et al. 20121. Chemical abundances 


can also provide clues on the relative importance of addi¬ 
tional evolutionary processes in galaxies, including stellar 
migration, gas flows within the disc and large scale galactic 
fountains. 

Current large integral field spectroscopy (IFS) galaxy 


surveys (eg. CALIFA, Sanchez et al. 2012 SAMI, Croom 
et al.|[20T^ MaNGA, Bundy et al.||2015 l offer great poten¬ 


tial for extending the study of resolved chemical abundances 
to much larger galaxy samples. Using optical IFS data, the 
chemical abundance of the interstellar medium (ISM) can 
be derived by analysing gas emission lines, while insight into 
the metallicity of the stellar population can be obtained us¬ 
ing Lick indices, or applying spectral decomposition tech¬ 
niques (Sanchez-Blazquez et al.]|2014 Cid Fernandes et al. 


2011). The CALIFA survey has recently delivered important 


insights on the distribution of metals within star-forming 
galaxies, providing the strongest evidence so far for a uni¬ 
versal gas phase abundance gradient ( [Sanchez et al.|[2014[ ) 
out to 2 effective radii (Re), as already suggested in pre¬ 
vious work (Vila-Costas & Edmunds 1992 Bresolin et al. 


|2009| [Bresolin, Kennicutt fc Ryan-Weber||20i(^ 

However, observations of resolved gas metallicity alone 
cannot be uniquely related to specific chemical evolution¬ 
ary models and, unless measurements of the gas and stellar 
content are also available, the information provided by the 
gas metallicity is highly degenerate between inflows and out¬ 
flows. 

This work presents one of the hrst attempts to com¬ 
bine chemical abundances from IFS with gas mass measure¬ 
ments from millimetre/radio observations for a single galaxy 
(NGC 628) on the same spatial scales, along with stellar 
mass surface density and stellar metallicities. We are there¬ 
fore able to constrain the metal budget both as a function of 
radius, and for the galaxy as a whole, out to 3 Ro. Since we 
have access to information on radial gradients of the relevant 
physical quantitates, we obtain an overall metal budget free 
from aperture effects, which affected previous work based on 
metallicities estimated from observations of galaxy central 


regions (for example from the SDSS galaxy sample, York 
[et al.|2000[ [Abazajian et al.|200^ . 

Empirical studies of the metal budget in the z = 0 Uni¬ 
verse, presented by several authors, agree in concluding that 
the total mass of metals produced by galaxies cannot be ac¬ 
counted for by the metals observed in stars and the ISM 


(|Pettini et al.[ 

1999 |Ferrara, Scannapieco & Bergeron|2005 

Bouche et al.| 

2007| Peeples & Shankar||2011| [Zahid et al.| 


2012 


Peeples et al. 20141. Large amounts of metals (between 


35% and 90%) must be residing in the intergalactic medium, 
but also in the diffuse halo gas and in the circumgalactic 
medium (CGM). Although recent observational efforts have 
led to a better characterisation of the metallicity of the gas 
in the hot halo phase | Tumlinson et al.|2011 20131, a robust 
metal budget for these phases is still missing. 

Cosmological hydrodynamical simulations which in¬ 
clude chemical evolution in a self-consistent way generally 
agree with the empirical studies above, and highlight that 
the fraction of missing metals is strongly dependant of the 



Figure 1. Optical image of NGC 628 (r' band from SDSS) with 
superimposed contours of Hi (from the THINGS survey, in blue) 
and CO surface brightness (from the HERACLES survey, in red). 
The red contours correspond to a molecular hydrogen surface 
mass density of 10 M 0 pc“^, while the blue Hi contour corre¬ 
spond to a neutral hydrogen surface mass density of 12 M©pc“^. 


sub-grid feedback prescription (Wiersma, Schaye & Theuns 
[2011[ [Pilkington et al(][2012[ ). 


The aim of this work is to investigate in detail how met¬ 
als are lost from the galaxy, hence the severity of the ‘missing 
metal problem’. This information also provides tight con¬ 
straints on the outflow loading factor (i.e. the ratio between 
SFR and outflow rate) averaged over the galaxy lifetime. 
The results are expected to provide simulators with new 
constraints on the net effect of gas flows in disc galaxies on 
resolved scales. 

The paper is structured as follows. In Sec. we give a 
brief overview of the observational data used, in Sec. we 
summarise how we use the data to derive gas phase metal¬ 
licity, gas and stellar mass surface densities and their radial 
gradients. In Sec. we summarise the chemical evolution 
framework, while in Sec. we present the results on the 
metal budget and the modelling of the outflow loading fac¬ 
tor in NGC 628. In Sec. and we present the discussion 
and conclusions. 


2 OVERVIEW OE THE DATA 


NGC 628 is a nearby (D ~ 7.3 Mpc, z = 0.00219) nearly 
face-on (i « 7°) spiral galaxy (Fig. [^. Table summarises 
some of its basic properties. NGC 628 is a good example 
of an isolated galaxy (no encounter in the last Gyr; see 


Kamphuis & Briggs 19921, displaying typical grand-design 


structure. The following sections describe the observations 
of NGC 628 that have been used in this work. 
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Galaxy gas flows and metal budget 3 


Table 1. General properties of NGC 628 


Property Value 

Name NGC 628. M 74. UGC 1149 


RA 

Dec 

Type 

Distance (adopted) 
i 

R24 

Re 

scale 


00'*59’"50.1^ 

-07°34'41.0" 

SAc 

7.3 Mpc 

70 

4.88' = 10.3 kpc 
67" = 2.4 kpc 
36 pc arcsec“^ 


2.1 Integral field spectroscopy 


In this study we make use of IFS observations of NGC 628 


taken as part of the PINGS survey (Rosales-Ortega et al. 


2010 Sanchez et al.|[?011 Rosales-Ortega et al.lpOlT l. The 

survey was carried out using the 3.5-m telescope at the Calar 
Alto observatory, with the Potsdam Multi-Aperture Spec¬ 
trograph in PPAK mode ( Kelz et al.|2006 1. The PPAK fibre 
bundle corresponds to a field of view of 75' x 65" and consists 
of 331 science fibres of 2.7" diameter packed in an hexagonal 
pattern, leading to a filling fraction of 65%. 

The data covers the spectral range between 3700 A - 
7100 A, with a spectral resolution of FWHM ~ SA. Due 
to the large projected size of the galaxy (10.5' x 9.5') the 
data consists of a mosaic of observations taken on 6 different 
nights over a period of 3 years. The hnal size of the mosaic 
is of 6' X 7'. We estimate the PSF to have a FWHM of 
approximately 3". 


2.2 HI 21 cm observations 


We used observations of the 21 cm line from THINGS (The 
Hi Nearby Galaxy Survey, Walter et al.||2008 l. This survey 
comprises 34 nearby galaxies observed with the VLA, at 
high spectral (< 5.3 km s“^) and spatial resolution. NGC 
628 was observed using the B, C and D configurations with 
a combined on-source integration time of ~ 10 hr. We use 
the publicly availably ‘robust’ weighted map (the ‘robust’ 
weighting scheme is used to achieve a sensitivity close to 
natural weighting while preserving a resolution close to uni¬ 
form weighting, see [Walter et aL]|2008| for details), which 
has a beam size of « 6 (Bmin = 5.57 ,Bmaj = 6.8"). The 
la noise per 2.6 km s“'^ channel is 0.66 mjy beam“'^, cor¬ 
responding to a sensitivity of « 0.5 Mq pc“^, sufficient to 
trace the atomic hydrogen in all regions where it constitutes 
the dominant component of the ISM. 


2.3 CO(2 -1) observations 

Observations of the CO(2 -1) transition we used to infer the 
molecular gas content. Maps are obtained from the publicly 
availablif] HERA CO-Extragalactic Survey (HERACLES, 


Leroy et al. 20091. HERACLES is a survey of 48 nearby 


galaxies using the HERA multi-pixel receiver on the IRAM 
30 m telescope, with 13" angular resolution and 2.6 km s“^ 


^ http://www.mpia.de/THINGS/Data.html 
^ http://www.mpia-hd.mpg.de/HERACLES/Data.html 


spectral resolution. The la sensitivity of the map is esti¬ 
mated to be ~ 3 Mq pc“^ (with the Galactic conversion 
factor discussed in Sect.3). 


2.4 Multi-Wavelength photometric data 

Photometry in different bands was collected to obtain a reli¬ 
able estimate of the stellar mass surface density by perform¬ 
ing pixel-by-pixel spectral energy distribution (SED) fitting. 
We collected observations from GALEX (Galaxy Evolution 
Explorer), 2MASS (Two Micron All Sky Survey) and the 
Spitzer Space Telescope IRAC camera. In the optical band, 
we make use of the photometry from SDSS Data Release 7. 
Table summarises the basic properties of the dataset used 
and the relevant references. 


3 DATA ANALYSIS 

In this section we present the analysis performed on the 
reduced data to obtain gas phase metallicities, stellar mass 
surface density and gas mass surface density profiles. We also 
briefly discuss the adopted choice of gas phase metallicity 
calibrators (Sec. |3.2[ ) and the radial gradients of the derived 
quantities (Sec. |3.6[ ). Throughout the section we highlight 
potential sources of systematic uncertainty. 


3.1 Emission line flnxes 


We extract emission line fluxes in each spaxel from the IFS 
datacube following a very similar procedure to |Belfiore et al.| 
(20151. Here we give a brief overview of the main steps in¬ 
volved. 


(i) For each spaxel, we fit a linear combination of single 
stellar population (SSP) templates, after correction for sys¬ 
temic velocity and instrumental dispersion, using Penalised 
Pixel Fitting (Cappellari & Emsellem 20041. In this work 
we used a grid of 29 SSP templates, generated by using 


Maraston fc Stromback (20111 models based on the empiri¬ 
cal STELIB spectral library ( Le Borgne et al.|2003 l. A spec¬ 
tral region of 600 km s“^ is masked around each emission 
line we wish to ht. Strong sky lines are also masked. We do 
not attempt to extract stellar population parameters from 
the SSP hts. Such a study has recently been carried out by 


Sanchez-Blazquez et al. (20141, and in Sec. 5.1 we make use 


of the stellar metallicity gradient derived in their work. 

(ii) For each spaxel, we subtract the stellar pop¬ 
ulation fit from the observed spectrum to obtain 
a pure emission line spectrum. This is then fitted 
with a set of Gaussian functions, one per emission 
line (using least-squares minimization). The doublets 
[OHIJAA4959, 5007 and [NIIJAA6548, 6584 are set to have 


the same velocity and velocity dispersion and the ratio of 
their intensities is fixed to the theoretical one. We calculate 
emission line fluxes by integrating the flux under the fitted 
Gaussians. 

(iii) We apply a signal to noise cut to the emission line 
maps generated above, imposing a S/N larger than 5 in 
Ha and H/3. This threshold is intentionally high to exclude 
regions of diffuse emission (mainly inter-arm regions) where 
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Table 2. Multi-wavelength observations of NGC 628 used in this work to supplement the IFS mosaic from|Rosales-Ortega et al.|l|2010ll. 


Telescope 

F liter/Receiver 

Agff [A] 

Bandwidth [A] 

PSF FWHM ["] 

Reference 


GALEX 

FUV 

1516 

268 

4.3 

Morrissey et al. 

(2007 



NUV 

2267 

732 

5.3 




SDSS 


u' 

3540 

599 

1.49 

g' 

4770 

1379 

1.34 

r' 

6222 

1382 

1.16 

i' 

7632 

1535 

1.03 

z' 

9060 

1370 

1.13 


Abazajian et al. (2009 I 


2MASS 

J 

H 

Kg 

1.235 [^lm\ 

1.662 [/.tm] 

2.159 

0.162 [/Am] 

0.251 [/Am] 

0.262 [/Am] 

2.5 

2.5 

2.5 

Skrutskie et al. (2006 1 

Spitzer 

IRACl 

IRAC2 

3.550 [/Am] 

4.493 [/Am] 

0.75 [pim] 

1.01 [pLm] 

1.66 

1.72 

Fazio et al. (2004 \ 

' - n —’ *-* 

IRAM 30m CO (2-1) 

HERA 

u = 230.54 [GHz] 


13 

Leroy et al. (2009 \ 


VLA Hi 


V = 1420.40 [MHz] 


5.57 X 6.80 


Walter et al. I 2008 1 


the Ha emission might not be ascribed to bona-hde Hll re¬ 
gions. We also discard regions where both the [Oil] and the 
[OIII] line are undetected. 

(iv) We calculate the reddening from the Balmer decre¬ 


of the most reliable, at least in the low metallicity regime 


ment, by using the Ha/H/3 ratio and a Calzetti (20011 at¬ 
tenuation curve with Ry = 4.05. The theoretical value for 
the Balmer line ratio is taken from lOsterbrock fc Ferlandl 
(20061, assuming case B recombination (Ha/H/I = 2.87). We 


note that the use of extinction curve of|Cardelli, Clayton fc| 
Mathis (19891 (or the modification by C)’Donnell|1994 1 with 


Ry = 3.1 yields very similar results for the 3600 A to 7000 
A wavelength range considered in this work. 


3.2 Gas phase metallicity 

Reliably measuring the gas metallicity from emission lines 
ratios remains a difficult problem in observational astro¬ 
physics. Different line ratios (‘diagnostics’) and methods 
(‘calibrations’) have been developed for the task. However 
it is well-known that different calibrations, even when based 
on the same diagnostics, can give results differing by up to 
0.6 dex ( Kewley fc Ellison||2008 Lopez-Sanchez et al.|2012 


[Pena Guerrero, Peimbert fc Peimbert|2012| . It is beyond the 
scope of this work to resolve the abundance scale problem 


(but see recent advances towards a possible solution in Do- 
|pita et al.|2013[[Perez-Montero|2014[|Blanc et al.|2015[ ). Fix¬ 

ing the gas phase metallicity scale is, however, a key factor 
in determining the total metal budget. We therefore discuss 
some of the assumptions going into widely used metallicity 
calibrators and possible shortcomings of each. 

The metallicity (i.e. the oxygen abundance 12 -|- 
log(0/H), where O and H are number densities), can be 
calculated knowing the electron temperature, which can be 
derived from the ratio of the oxygen auroral line [OHI]A4363 
and another oxygen line, like [OIII] A5007. This method (re¬ 
ferred to as the Te method^ is generally considered one 

^ Or sometimes as the ‘direct method’. However, given the well- 
known systematic uncertainties in Tg-measured metallicities, we 
feel that this designation is inappropriate 


(Pagel et al. 1992 Garnett 1992 Izotov et al. 20061. At 


high metallicities, however, the electron temperature de¬ 
creases and the auroral line is generally undetected. More¬ 
over, temperature fluctuations make the Tg method unre¬ 
liable at high metallicities | Stasins^|2005 Bresolin]|2007 1. 
Due to the strong temperature dependance of the emissivity 
of the auroral lines, regions with higher than average tem¬ 
perature may dominate the line emission, hence biasing the 
metallicity measurements towards low values. Finally, there 
is evidence that the distribution of electron energies in Hll 


regions might not follow a Boltzmann distribution (Binette 


et al. 2012 Nicholls et al. 20131, hence making the defini¬ 


tion of temperature arguable, while a ^-distribution may be 
more appropriate (Dopita et al.|20T3|. 


Recombination lines can also be used to estimate the 
metallicity. These are almost insensitive to temperature, 
which makes them much less dependent on temperature 
fluctuations and on the specific distribution of electron en¬ 
ergies. However, these lines are approximately 10"* times 
fainter than H/3, which makes their detection possible only 
for bright sources on 8-meter class telescopes. Metallicities 
estimated from the recombination lines are systematically 
offset, by about 0.2 dex, compared to those measured with 


the Tg method ( 

Blanc et al. 

2015 

1. Some authors ( 

Peim- 

bert & Peimbert 

|2010| |Pena Guerrero, Peimbert & Peim- 

bert 2012 

1 attribute this discrepancy to the presence of small 


temperature fluctuations within Hll regions, which (as dis¬ 
cussed above) would bias the observed auroral line fluxes 
towards the highest temperature, lower abundance regions. 

Finally, several authors have calibrated 
strong nebular line ratios (making use of 
[OHI]AA4959,5007, ]OH]A3727, [NH]A6584, etc) as abun¬ 
dance diagnostics. While the intensity of these lines depends 
on several factors besides the metal abundance (ionization 
parameter, density, N/O ratio), calibrations take advantage 
of the fact that in the Hll regions observed in the local 
Universe relations exist between some of these parameters. 
Strong line diagnostics can be calibrated against Tg-based 
metallicity measurements (Pilyugin, Vflchez & Thuan 
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2010 Pilyugin & Thuan 2005 Pettini & Pagel 20041, 


photoionisation models (Denicolo, Terlevich & Terlevich 


2002 Tremonti et al. 2004 Kobulnicky & Kewley 1 2004 1, 
or a mixture of both (|Nagao, Maiolino fc Marconi| 


Maiolino et al. 20081. Recently, Blanc et al. (20151 has 


2006 


shown that some (but not all) photoionisation model grids 
are able to reproduce the metallicities measured using 
recombination lines. 

To take into account the systematic uncertainty intro¬ 
duced by the choice of metallicity calibration, in this work 
we use two independently derived calibrations. 


(i) The calibration from Maiolino et al. 

(2008|) (M08), 

which is based on photoionisation models 

Kewley & Do- 


pita 2002 \ in the high metallicity regime and anchored to 


Te-based metallicity measurements in the low metallicity 
regime. While M08 provides calibrations for various strong 
line diagnostics, here we use only the calibration based on 
the R23 parameter, 

R 23 = ([OII]A3727 + [OIII]A4959 -b [OIII]A5007)/H,d. 

In the context of this work we have re-calculated the 
parametrisation of the R23-metallicity relation using the 
same procedure in M08, but using a 5th order polynomial 
fit (M08 used a fourth order polynomial) to provide a bet¬ 
ter fit to the high-metallicity end of the relation between 
12-|-log(0/H) and R23. Overall the M08 calibration gives 
similar abundances to those obtained using the calibration 
used by [Tremonti et ah] (|2004[) and the calibration of R23 
presented by |Kobulnicky fc Kewley (20041. Since the R23 
parameter is double-valued, we break the degeneracy by us¬ 
ing the [Nil]/[Oil] ratio, concluding that all the regions in 
NGC 628 belong to the upper branch of R23. We note that 
the use of a calibration the recursively solves for both metal¬ 
licity and the ionisation parameter ( [Kobulnicky fc Kewley] 
2004) does not substantially alter the derived metallicities. 


(ii) The calibration from Pettini & Pagel (2004) (PP04), 
based mostly on Te measurements, using the 03N2 index: 


03N2 = log 


[OIII]A5007/H,a 

[NII]A6584/Ha 


This calibration has the disadvantage of depending explicitly 
on the the nitrogen abundance, which does not necessarily 


scale in a simple way with the oxygen abundance (Perez- 


Montero fc Contini|2009||Perez-Montero|2014[|Belfiore et~ 


2015), since nitrogen has both a primary and a secondary 


nucleosynthetic origin and is released into the ISM on longer 
timescales than oxygen. 

It is also worth noting that we are (perhaps unduly) ex¬ 
tending the PP04 calibration to super-solar metallicities. In 
the original work of Pettini & Pagel (2004), only two Hll re¬ 


gions with super-solar metallicity and Te measurements are 
presented, while the high-metallicity end of the calibration 
is fixed by considering a small number of regions with metal¬ 
licities calculated using photoionisation models. Even in the 
most recent calibration of the 03N2 diagnostic, presented 
by Marino et al. (2013), only 5 Hll regions with super-solar 


metallicity are included, which makes any attempt to ‘di¬ 
rectly’ calibrate the relation between 03N2 and metallicity 
in the super-solar regime subject to large extrapolation un¬ 
certainties. 

We also note that the two adopted calibrations make 


different assumptions regarding dust depletion of oxygen. 
In the M08 calibration, an oxygen depletion factor of 0.22 


dex is assumed (Kewley & Dopita 2002), while Te-based 


calibrations do not correct for the amount of oxygen de¬ 


pleted onto dust. We follow Peimbert & Peimbert (2010) 


and Pena Guerrero, Peimbert & Peimbert (2012) and mod¬ 


ify the metallicities obtained from both calibrations assum¬ 
ing a constant oxygen depletion factor of 0.1 dex. For the 
rest of this work we will use gas phase metallicity to refer 
to the overall oxygen abundance of the ISM (gas and dust), 
with the understanding that the dust is taken into account 
through a constant depletion factor. Altering the dust de¬ 
pletion by ±0.1 dex does not substantially change any of our 
conclusions. 

Since M08 and PP04 are representative of the two 
main classes of metallicity calibrators (Te-based and 
photoionisation-model-based), alternative calibrations will 
not lead to metallicity estimates significantly higher than 
M08 or significantly lower than PP04. A thorough analy¬ 
sis of the gas phase metallicity in NGC 628 was presented 
in [Rosales-Ortega et al.j ( |2011| ) , who used a similar set of 
calibrators (as noted above, the Kobulnicky fc Kewley|2004 


calibration is roughly equivalent to the adopted [Maiolino 


[et al.[[^08P and obtained similar gradients and zero-points 
to the ones presented in this work. Interestingly, recent work 


from Croxall et al. (2013), who used far-IR lines to estimate 


metallicity in NGC 628, points towards oxygen abundances 
within the range bracketed by the M08 and PP04 calibra¬ 


tions. Additionally Berg et al. (2015) obtained direct T, 


measurements of metallicity for 45 Hll regions in NGC 628, 
which agree on average, though with large scatter, with the 
metallicities derived with the PP04 calibration. The reason 
for the large scatter in the Tg-derived metallicities is not 


clear, but according to Berg et al. (2015) might be due to 


systematics in their electron temperature measurements. 

Since metallicity diagnostics are only calibrated for star 
forming regions (classical Hll regions), regions of galaxies 
where other types of ionisation (shocks, active galactic nu¬ 
clei, evolved stars, etc) are dominant should be excluded. 
We use the standard ionisation diagnostic diagram (BPT 
diagram Baldwin, Phillips fc Terlevich] 198 1] [Veilleux fc Os-[ 
[terbrock 1987 Kewley et al. 2001] Kauffmann et al. 2003) 
to identify galactic regions as star-forming and discard re¬ 
gions that lie above the [Kewley et a~ (2001) demarcation 
line in the [OIII]/H/3 vs [SII]/Hq diagnostic diagram. All 
of the discarded spaxels are found in the inter-arm regions 
and present LINER-like ionisation, which might be due to 
evolved hot stars dominating the ionisation budget in re¬ 


gions where no recent star formation has taken place (Yan 
& Blanton 2012 Belfiore et al. [2015 ). We note that such 


regions, not complying with the Hll classification, are only 
5% of the total, and their exclusion does not affect our final 
results. 

The resulting maps of gas phase metallicity for the M08 
and PP04 calibrations are shown in Fig. For reference, 
the photospheric solar oxygen abundance, as derived by[As-[ 


plund et al. (2009), is 12±log(0/H) = 8.69, or equivalently 


^©(O) = 0.00585. Note that we are referring to Z as the 
oxygen abundance by mass. The total metal content (frac¬ 
tion by mass of elements heavier than helium) in the Sun is 
0.0142. 
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12 + log(0/H) M08 



12 + log(0/H) PP04 



Figure 2. Maps of the gas metallicity in NGC 628 using the |Maiolino et al.| | |200^ (M08) calibration based on R23 (left) and the |Pettini| 
I&: Pagei| l |2004| l (PP04) calibration based on 03N2 (right). Only spaxels that lie below the |Kewley et aT] | |200ll l demarcation line in the 
[OIII]/H/3 versus [SII]/Ha BPT diagram are included. 
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Figure 3. Left: Stellar mass surface density map, created from the spatially resolved SED fitting. Foreground stars have been masked. 
Right: Radial dependance of various physical quantities across the disk of NGC 628 (in radial bins; error bars represent the scatter 
in each radial bin). Atomic gas surface density (blue) shows a dip in the central region, while the molecular gas surface density (red) 
is peaked at the centre and quickly decreases outwards. The horizontal red dashed line corresponds to the sensitivity limit of the CO 
observations. The total gas content (green) is nearly fiat across the whole disk up to 3 Re (1 Re = 2.4kpc). A clear radial gradient is 
observed in both stellar mass surface density (black, short dashes) and gas metallicity (violet). The latter is shown both for the M08 
(solid violet) and PP04 (dot-dashed violet) metallicity calibrations. The stellar metallicity gradient from |Sanchez-Blazquez et al.|p014[ | 
is shown in orange. Note that the metallicities of gas and stars are plotted with respect to the alternative y-axis on the right in units of 
12-hlog(0/H). 
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Galaxy gas flows and metal budget 7 


3.3 Stellar mass surface density 

Stellar mass surface densities are estimated by performing 
SED fitting in each pixel on the multi-wavelength photo¬ 
metric data from UV to near IR, using the software package 


Worthey et al. 19941, which acts in the direction of con- 


CIGALE (Code Investigating GAlaxy Emission, Burgarella, 
Buat fc Iglesias-Paramo|2005 Noll et al.|2009| . In detail, we 
followed the procedure outlined below (similar to Boquien 
et al.|2012| ). 


(i) The different observations were smoothed to the low¬ 
est resolution of the photometric dataset (GALEX NUV), 
using the kernels from [Aniano et al.| ( |2011[ ) . 

(ii) The maps were then regridded on the same pixel scale 
as the IFS mosaic. 

(iii) For each pixel, the flux and the flux error in each 
band was extracted. Errors are computed by taking into ac¬ 
count the uncertainty in absolute calibration of the different 
instruments. Fluxes and errors are used as input for the 
software package CIGALE. 

(iv) CIGALE creates FUV-to-FIR SEDs consisting of dust- 
attenuated complex stellar population models. In this anal¬ 
ysis we used [Maraston] ( |2005[ ) templates (with a Kroupa 
IMF) to build a grid of stellar models with a range in metal- 
licity (from 0.5 to 2 solar). The star formation history is 
parametrised by an exponentially decreasing SFR, and an 
old stellar population with and a grid of younger stellar pop¬ 
ulation bursts. Extinction is fitted using a Calzetti-like at¬ 
tenuation curve. The optimal galaxy parameters are derived 
by CIGALE using a Bayesian-like analysis. A map of stellar 
mass surface density is then created (Fig. left panel) 

The formal errors computed by CIGALE (which take into 
account the uncertainty in the data, but not the systematics 
of the method) are of the order of 0.1 dex. As a further check, 
we compared the stellar mass estimated by CIGALE with the 
stellar mass estimate obtained using the procedure in |Bell| 
et al. (20031 and using the u-r colour to estimate the mass- 


to-light ratio. The stellar mass densities obtained with the 
Bell formulation are consistent with the ones obtained with 
CIGALE, with a small systematic difference of less than 0.1 
dex at very large radii (R > 8 kpc). The total mass of NGC 
628 out to 3 Re computed by CIGALE is log(M*/M0) = 9.8. 


3.4 Stellar metallicity 

To take into account the metals locked in stars we make 


use of the recent study of Sanchez-Blazquez et al. (20141, 


who present the mass-weighted stellar metallicity gradient 
for NGC 628 by performing full spectral fitting using the 
software package steckmap on the PINGS IFS data. De¬ 
spite the increased level of sophistication of the spectral fit- 


ting algorithms in recent years 1 Heavens, Jimenez & Lahav 

20001 [Ranter, Heavens & Jimenez||2004| 

Fernandes, Sodr & 

Gomesj 2005| Gallazzi et al.j 20061 [Sane 

lez-Blazquez et al. 

2011| Gid Fernandes et al.|2014||Wilkinson et al.|2015 [Me- 

Dermid et al.||2015 

1, stellar metallicity remains the hardest 


quantity to estimate in stellar population studies. The diffi¬ 
culty is due to the non-linear effect of a limited set of metal¬ 


fusing old, metal poor systems with young metal rich ones. 
While the ability of spectral fitting codes to reliably derive 
input parameters increases as a function of S/N, it has been 
shown (see for example Cid Fernandes et al.|2014| Wilkinson 


et al. 20151 that even for high S/N spectra, stellar metal¬ 


licity calculated with different sets of simple stellar popula¬ 
tions can differ by more than 0.2 dex. It is beyond the scope 
of this work to perform a new full spectral fitting study of 
the PINGS data for NGC 628 and further evaluate these 
systematic uncertainties. We note, however, that our simple 
approach of propagating of the errors in the stellar metal- 


licitles presented in Sanchez-Blazquez et al. (20141 (median 
error ~ 0.15 dex) might represent an underestimate of the 
real uncertainty. 

Taking into account the necessary correction factor due 


to the fact that Sanchez-Blazquez et al. (20141 adopt a dif¬ 


ferent distance for NGC 628, we can obtain a stellar metallic¬ 
ity gradient out to 5 kpc. We put the gas phase and stellar 
metallicity on the same abundance scale by using the so¬ 


lar metallicity of Asplund et al. (20091 and assuming solar 


abundance ratios. We do not attempt to correct for possi¬ 
ble non-solar abundance ratios (e.g. a-enhancement) in the 
bulge, since in any case the correction will only affect our 
innermost radial bin and will not have a significant effect on 
the subsequent analysis of integrated properties. 


3.5 Gas mass surface density 

The Hi mass can be calculated directly from the Hi surface 
brightness, assuming that Hi is everywhere optically thin. 

To convert from GO luminosity to H2 mass we adopt a 
constant CO to H2 (aco) conversion factor given by 


aco = 4.35 


M0 pc 


K km s-i 

This value of the conversion factor is appropriate for the 


Milky Way (Solomon et al. 1987 Strong & Mattox 1996 


|Abdo et al.||2010[ [Bolatto, Wolfire fc Leroy||2013[ ). We ex¬ 
plored the effect of using a metallicity-dependant conversion 


factor (Schruba et al. 20111 and the conversion factor de¬ 
rived by [Blanc et al.| ( |2013[ ) by inverting the star-formation 
law in NGC 628, finding that none of our conclusions are 
substantially modified. As the conversion factors are cali¬ 
brated for the CO (1-0) line, while the HERACLES obser¬ 
vations map the CO (2-1) transition, we assume a constant 


ratio of 0.7 between the two CO luminosities (Leroy et al. 

[M2| |. 

To calculate the total gas mass we multiply the Hi com¬ 
ponent by 1.36 to account for helium and add the molecular 
gas mass computed using the conversion factor above, which 
already takes helium into account. 


3.6 Radial gradients 

It is well known that spiral galaxies exhibit radial gradients 
in their physical properties, including stellar mass surface 


sented) and the well-known age-metalllcity degeneracy (e.g. 


licity values generally available in spectral libraries (since 

density and metallicity ( 

Bell & de Jong 

20011 

Moustakas 

most stellar libraries rely on observations in the solar neigh- 

et al. 1120101 [Gonzalez Delgado et al.||2013[ [Sanchez et al. 

bourhood super-solar metallicity stars are not well repre- 

20141. 


In Fig.j^we present the radial gradients for the different 
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8 F. Belfiore et al. 


physical quantities that we have derived in this section for 
NGC 628. We observe an exponential profile in stellar mass 
surface density and gas phase metallicity (black dashed and 
violet lines respectively). 

We note here that the PP04 gas metallicity gradient 
presents a flattening in the inner 2.5 kpc (« 1 Re, violet dot- 
dashed line), already reported in previous work ( [Rosales^ 
Ortega et al. 2011| and consistent with the flattening ob¬ 


served in a subsample of galaxies from the CALIFA survey 
| |Sanchez et al.|2014p . With the current data it is not possible 
to assess whether this feature corresponds to a real abun¬ 
dance drop in the central region or if it is an artefact of the 
metallicity calibration. In particular, we note that such flat¬ 
tening observed when using the 03N2 calibration is often as¬ 
sociated with a central decrease of the N/0 abundance ratio 
(which can be traced by the [Nil]/[Oil] flux ratio), and this 
is also the case of NGC 628. Therefore the metallicity flat¬ 
tening may actually reflect a relative abundance variation af¬ 
fecting the metallicity diagnostic used in PP04, rather than 
a global metallicity flattening. However, we note that from 
the theoretical standpoint, bars and other non-axisymmetric 
perturbations can trigger large scale gas flows and induce a 
flattening of the metallicity gradient (Roskar et al.|2008a|b 


Cavichia et al. 20131, together with a metallicity ‘hump’ 


at the interface. While the gas phase abundance does not 
show evidence for this hump, a sharp increase in the stellar 
metallicity is reported at 2.5 kpc by |Sanchez-Blazquez et ah] 
(20141 (orange line in Fig. [^. They interpret this feature 
as the effect of an oval distortion (possibly a signature of 
a dissolving bar). The presence of a circum-nuclear ring of 
enhanced star formation seems to confirm this hypothesis. 

The Hi disc is much more extended than the FoV of the 
IFS data and presents clumpy substructure (seen in Fig.j^. 
The Hi surface density profile is characterised by a central 
dip and and an increase towards larger radii. The maximum 
in Hi surface density coincides roughly with the edge of the 
IFS field of view (« 3 Re), hence in our subsequent analysis, 
which is limited to the radial extent of the IFS observations, 
we neglect about 40% of the mass of the Hi disc. We will 
discuss the importance of this extended Hi disc in Sec. 

The molecular gas profile is centrally concentrated and 
its surface density decreases with radius. The sensitivity 
limit of the HERACLES data is shown as a horizontal dotted 
line in Fig. The combination of the atomic and molecu¬ 
lar gas profile generates a total gas surface density profile 
(green line in Fig. which is remarkably constant over the 
whole field of view of the IFS data. This implies that the 
gas fraction (fgas = Mgaa/(M* -|- Mgas)) increases outwards 
roughly log-linearly with radius. 


4 THE CHEMICAL EVOLUTION 
FRAMEWORK 

In this section we present the analytical framework we will 
be using to study the metallicity evolution in resolved re¬ 
gions of NGC 628. 

We make use of simple relations obtained by invoking 
the instantaneous recycling approximation and perfect in¬ 
stantaneous mixing. This framework has been proven suc¬ 
cessful in modelling the chemical evolution of oxygen (Zahid 


|et al.|[^012| [Dayal, Ferrara fc Dunlop||2013[ [Peeples et al. 


2014 Peng fc Maiolin^|2014 1, which is predominantly pro¬ 
duced by short-lived massive stars dying as core-collapse su¬ 
pernovae. We note in passing that the chemical abundances 
of other common elements (like N or Fe) require a more de¬ 
tailed model, taking into account stellar lifetimes, the Type 
la supernovae delay time distribution and variation of the 
nucleosynthetic yields with metallicity, and will not be at¬ 
tempted in this work. However, since oxygen is the most 
abundant metal by mass, it is a good tracer for the total 
metal content. 

Several studies in the literature have shown that, in 
order to interpret observed chemical abundances in galax¬ 
ies, it is necessary to devise a model taking gas inflows and 
outflows into account (often dubbed the ‘gas regulatory’ or 
‘bathtub’ models; |Matteuca]|1986||Gibson et al.|2003||LilIy| 


et al. 2013 Dekel fc Mandelker|2014[ Peng fc Maiolino|2014 1 . 

Within the instantaneous recycling approximation, de¬ 
noting the oxygen fraction (by mass) in the ISM as Zg, the 
gas mass as Mg, the star formation rate as SFR, the stellar 
mass as M*, the total mass in oxygen in the galaxy as Mz, 
the mass of oxygen in the ISM as Mzg, the outflow rate as Vli 
and the inflow rate as 4?, the galaxy’s chemical evolution in 
a bathtub model is described by the following constitutive 
equations 


^ = -F - (1 - R) SFR - 4-, 


dM„ 

dt 


= (1 - R) SFR, 


( 1 ) 


( 2 ) 


dMzg ^g) 


dt 

and 


dt 


= p SFR - ^Zg - (1 - R) Zg SFR, (3) 


dMz 

dt 


= p SFR - ^Zg 


(4) 


where p is the mass of newly synthesised oxygen per 
unit mass of gas converted into stars (we will refer to it 
as the yield per stellar generation) and R is the fraction of 
gas mass that is promptly returned to the ISM via stellar 
mass loss processes (referred to as the return fraction). We 
emphasise that the same relations apply also to the total 
abundance of metals (after replacing the appropriate value 
for the net nucleosynthetic yield), if the abundance ratios of 
different elements are assumed constant. 

Both the average yield and the return fraction are a 
function of the initial mass function (IMF) and in general 
also depend on metallicity and time. In our simple model 
neither the time nor the metallicity dependance of R and y 
will be further considered, since the oxygen yield is shown to 
be approximately independent of metallicity and the effect of 
stellar lifetimes on the return fraction can be considered neg¬ 
ligible for our purposes ([Thomas, Greggio &; Bender 1998 


Kobayashi et al.||2006[ [Zahid et al.||2012[ [Vincenzo et al.| 
2015 ^ 


If we denote as p(M) the net yield of oxygen for stars 
of mass M, then the net yield (p) per stellar generation is 
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Galaxy gas flows and metal budget 9 


given by 


P = 


IMF(M) dM 


/^3IMF(M) dM 


(5) 


where the IMF(M) is the initial mass function, Miow and 
Mup are the lower and upper mass cutoffs of the IMF, while 
Miog-iiv is the highest mass of the stars that are consid¬ 
ered ‘eternal’ or, equivalently, the lowest mass of the stars 
contributing to the chemical enrichment. Note that an alter¬ 
native dehnition of yield is given by the new mass of oxygen 
produced per unit mass of long-lived stars l[Searle & Sar- 


(m this case generally denoted as y). These two 


definitions are trivially related by the return fraction with 
y — p/(l — /?)|^It is important to stress that the yield per 
stellar generation is strongly dependent on the IMF, and 
not only on the ‘stellar’ nucleosynthetic yield p(M). In fact, 
the average yield can vary by a factor larger than three for 
the same choice of stellar yields if the IMF is changed from 


Kroupa, Tout & Gilmore (19931 to Chabrier {20031 (Vin- 
et al.|2015D. 


In this work, when a numerical value for the average 
yield is required, we make use of the results in |Vincenz^ 
et al. (20151, based on the compilation of stellar yields in 
Romano et al.| ( |2010[ ), which have been shown to success¬ 
fully reproduce the oxygen abundance in the Milky Way. 


We adopt a KTG IMF (Kroupa, Tout & Gilmore 19931, 


Miong-iiv = 1 Mq and therefore get a fiducial average yield 


of p = po = 0.007 for oxygen and R = 0.30 (Vincenzo et al. 


2015). The same results can be applied to the total content 


of metals, by replacing the numerical value of the average 
yield with p = pz(tot) = 0.013. 

It can be easily shown that in the case of a ‘closed box’ 
evolution, i.e. no inflow and no outflow (<0 = 0 and ^' = 0), 
the gas metallicity is given by 


Z, = 


1 - R 


In (fgas) = yin (fgas), 


( 6 ) 


where fgas = Mg/(Mg-|-M*) is the gas fraction. Some au¬ 
thors find convenient to define the so-called ‘effective yield’ 
as 


Jeff = 


In(fg-al) 


(7) 


Whenever yefi 7^y = p/(l — R) Eq[^ implies that the 
system has not evolved as a closed box. Fig. shows the 
radial gradient of the effective yield in NGC 628, clearly 
demonstrating that increases with radius. Importantly, 
the observed effective yield is always lower than y = p/(l — 
R) = 0.010, implying that the region of the disc of NGC 628 
sampled by us has not been evolving as a closed box. 

However, the information provided by the effective yield 
is highly degenerate. More specifically, both enriched gas 
outflows or pristine (or low metallicity) gas inflows will lead 
to the observed effective yield to be lower than the closed 
box value |Edmunds|1990|). 


^ Note that |Peng &; Maiolino| (|2014|| does not follow this tradi¬ 
tional notation. The quantity denoted as y in |Peng fc Maiolin^ 
(2014 1 is the same quantity which we denote as p in this work. 



R [kpc] 

Figure 4. The radial variation of the effective yield in NGC 628. 
The value expected for the closed box case, y/(l — R) = 0.010, is 
indicated as a black dotted line. The red circles and blue triangles 
represent respectively the results obtained with the |Maiolino et ah] 
(|2008[l and [Pettini fc Pagel| |[2C)04[) calibrations. 


5 RESULTS 

5.1 A spatially resolved oxygen mass budget 

The aim of this section is to compute the total mass of oxy¬ 
gen present in NGC 628 as a function of radius and compare 
it with the total mass of oxygen expected to have been pro¬ 
duced by its stellar component. 

This comparison provides a direct (model-independent) 
measurement of the mass of metals lost by each region in 
the galaxy via metal-enriched outflows. Moreover, the re¬ 
sults obtained in this section can be directly compared with 
predictions from cosmological simulations. 

Following Eq. the total mass of oxygen produced is 
given by 


=y M*. 


( 8 ) 


This corresponds to the total mass of oxygen actually 
present in the galaxy in the case of no outflows (Eq. [^with 
fk = 0). The total observed mass of oxygen is given by the 
sum of the mass of oxygen in the ISM (Mzg = ZgMg) and 
the mass of oxygen locked in stars (Mz* = Z*M*), i.e. 

Mz(obs,tot) = Mzg -l- Mz* = ZgMg -|- Z*M* . (9) 


Both simulations (Marinacci et al. 2014) and observa¬ 


tions ( Shull, Smith fc Danforth]2012 Tumlinson et al.|2013 ) 
argue that a large quantity of baryons and metals might 
reside in the hot haloes around galaxies. However, for the 
purposes of chemical evolution this gas is irrelevant since we 
assume it does not form stars, at least as long as it remains 
in the hot phase. Therefore in this work we consider this 
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Figure 5. Left: Oxygen mass per unit area in different galactic components compared with the oxygen mass expected to have been 
produced by the observed stars, in different radial bins. The orange line shows the oxygen mass locked in stars, by using the stellar 
metallicity gradient from [Sanchez-BIazquez et al.| ( [2014[ l. The violet solid line shows the oxygen mass in the ISM obtained by using the 
M08 calibration, while the violet-dashed line shows the oxygen mass in the ISM obtained by using the PP04 calibration. In red we show 
the total observed mass of oxygen (ISM + stars). The blue line shows the mass of oxygen produced by the stars at the same location 
assuming our fiducial average yield. Top Right: Mass of oxygen in each radial bin as a percentage of the mass of oxygen produced in the 
same radial bin. The colour-coding is the same as in the left panel. Bottom Right: The cumulative mass of oxygen within radius R in 
the different galactic component, expressed as a percentage of the total mass of oxygen produced within the same radius R. The points 
shown for the last radial bin correspond to the integrated oxygen budget for NGC 628, showing that the galaxy as a whole has retained 
only ^ 50% of the metals it has produced. 


gaseous component as effectively external to the galaxy and 
the metals present in it as lost. 

In Fig. (left panel) we show the total observed mass 
of oxygen per unit area in stars (orange), in the ISM (vio¬ 
let), in ISM-bstars (red) and the mass of metals that must 
have been produced by the observed stellar mass using our 
adopted yield (using Eq.[^ blue line) as a function of radius. 
In the case of the ISM and ISM-|-stars lines, we show the re¬ 
sults obtained with both the M08 (solid) and PP04 (dashed) 
metallicity calibrations. The error bars represent the error 
in the mean in each annulus, added in quadrature to the me¬ 
dian of the intrinsic error in each of the physical quantities. 
The intrinsic error budget is dominated by the error on the 
stellar metallicity, followed by the uncertainty in metallicity 
and stellar mass (for both we get a median error of about 
0.1 dex). These errors do not attempt to incorporate system¬ 
atic effects, which would dominate the uncertainty budget 
for this study. Moreover we take the net yield (y) and the 
return fraction (R) to be known exactly. With these caveats. 
Fig. [5] demonstrates that using the fiducial yield at all radii 
more metals have been produced than can be accounted for 
considering both the metals in the ISM and in stars. 

In the discussion (Sec.[^ we argue that a different choice 


of IMF (adopting for example Salpeter or Chabrier IMF) 
only makes the metal deficit worse. We also note that the 
mass of metals in stars presents a sharp feature between 2.5 
and 3 kpc, which is a direct consequence of the sharp feature 
in the stellar metallicity gradient. Since it is present in only 
one radial bin, the effect of this feature on the overall metal 
budget of the galaxy is negligible and will not be discussed 
further. 

In Fig. 1^ top right panel, we show the mass of oxygen 
in stars (orange) and ISM (violet) as a function of radius as 
a percentage of the total mass of oxygen produced in that 
annulus. The total mass of oxygen in each annulus (ISM 
+ stars, red line) is also shown, and we plot the results 
using both the M08 (solid line) and PP04 (dot-dashed line) 
metallicity calibrations. 

Fig. Ill (top-right) shows that the fraction of metals lost 
is about 45%-55%, and approximatively constant within the 
uncertainties for R > 3 kpc. There is also some evidence that 
the fraction of metals lost in the central region is higher (only 
~30% of the metals have been retained). This feature may 
trace the effect of massive outflows during the formation 
of the bulge or the cumulative effect of nuclear activity in 
the past. Interestingly Fig. (top-right) also shows that the 
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Table 3. The total oxygen mass in different galactic components, 
integrating the radial gradients for NGC 628 out to R = 7.00 kpc. 


Metals 

Mass [lO"^ Mq] 

Mass [% Mz prod] 

Stars 

2.6 ±0.9 

25 ± 2% 

ISM (M08) 

3.2 ±0.9 

31 ± 3% 

ISM (PP04) 

2.3 ±0.7 

23 ± 2% 

Stars + ISM (M08) 

5.8 ± 1.7 

56 ± 5% 

Stars -1- ISM (PP04) 

4.9 ± 1.5 

47 ± 5% 


fraction of metals retained increases at R > 5 kpc. Gas phase 
metallicity measurements at even larger radii provided by 
the Berg et al.| (20151 dataset confirm and further extend 


this trend (see discussion in Sec. |6.2[ |. 

Fig. [^(lower-right) shows the cumulative version of the 
previous plot. For each radius R, the mass of oxygen in dif¬ 
ferent components within radius R is computed as a per¬ 
centage of the mass of metals produced within the consid¬ 
ered radius. This plot can be interpreted as the results that 
would be obtained for the integrated properties of NGC 628 
if all quantities (metallicities, gas masses etc) were measured 
using an aperture corresponding to a radius R on the sky. 

We note that cumulatively the mass of oxygen in the 
ISM is sub-dominant out to R « 5.5 kpc, at the latter ra¬ 
dius stars and ISM contribute roughly equally to the total 
oxygen budget. For easy reference, the masses of oxygen in 
the different components within 7.0 kpc are summarised in 
TablejS] Taking into account the uncertainty associated with 
the choice of metallicity calibration, we conclude that out to 
R = 7.0 kpc the bulge and the disc of NGC 628 as a whole 
have lost ~ 50% of the oxygen they have produced. 


5.2 Inferred net outflow loading factor 


As discussed above, in the absence of outflows the total 
amount of metals produced should be simply given by the 
stellar mass times the amount of metals produced per mass 
of long lived stars (Eqj^. To model the effect of outflows 
we follow the formalism laid out in Peng & Maiolino (20141 
and assume the outflow rate fk to be proportional to the star 
formation rate through a constant outflow loading factor A, 


4- = A SFR. 


( 10 ) 


The resulting total amount of metals is given by the equation 

Mz(tot) =yM,-xJ Zg(t) SFR(t) dt. (11) 

The last term is the ‘metal deficit’ and its evaluation 
requires knowledge of the evolution of the metallicity with 
the star formation rate. To make further progress we there¬ 
fore need to solve Eq. [TJ|^ A general solution as a function 
of time can be obtained assuming a star formation law. Fol¬ 
lowing previous work we use a simple linear relation of the 
form 


SFR = e M,j 


( 12 ) 


where e is often refereed to as the star formation efficiency 
or also as the inverse of the gas depletion time (l/ra). 

The final piece of information needed is the functional 
form of gas inflow rate il?. A number of authors have assumed 


a gas inflow rate proportional to the SFR (e.g. Matteucci & 



Figure 6. Metals-to-stellar mass ratio (Mz/M*) versus gas frac¬ 
tion (fgas). Symbols show the values observed in radial annuli in 
NGC628, where the colour-coding gives the galactocentric radial 
distance, as indicated in the colour bar. Solid symbols correspond 
to the mass of metals in which the ISM metallicity is calculated 
by using the M08 calibration, while open symbols use the PP04 
calibration. The red lines correspond to the prediction of simple 
chemical evolution models with different values for the average 
outflow loading factor A. Note that the model relations are inde¬ 
pendent of the inflow rate <I> and of the star formation efficiency 
£. 


Chiosi|1983||Erb|2008||Mannucci et al.|2009||Dayal, Ferrara| 


fc Dunlop|2013 Troncoso et al.|2014 Kudritzki et al.|2015 |. 

This assumption leads to simple analytical solutions of the 
resulting differential equations. Moreover it does not require 
consideration of the time variable, since all the quantities 
evolve in time like the SFR and hence the time dependance 
factors out. However, the assumption of an inflow rate pro¬ 
portional to the SFR has no physical grounds and has been 
used in previous work entirely for convenience. A much more 
realistic assumption is that of a constant, or slowly varying. 


inflow rate. As discussed in Peng & Maiolino (20141 this 


assumption is much closer to the expectation based on the 
inflow rate of baryons in dark halos inferred in analytical cal¬ 
culations and numerical simulations of structure formation 


( Faucher-Giguere, Keres fc Ma|2011 Dekel et al.|2013 1. 

Interestingly, analytical solutions of Equations [fjj^ can 
still be obtained with the assumption of constant inflow rate, 
although they do depend explicitly on time. Exact solutions 
for Mg and M* are given in |Peng fc Maiolino (20141. Further 
exact solutions for Zg and Mz/M* are presented in Appendix 
[A|of this work. Here, we exploit those solutions to derive the 
relation between Mz/M* and the gas fraction. We refer the 
interested reader to Appendix [^ for a detailed discussion of 
the assumptions, strengths and weaknesses of the different 
bathtub models and a justification of the choice made in this 
work. 
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We further show in appendix that the relation be¬ 
tween Mz/M* and fgaa is independent of both the inflow 
rate and of star formation efficiency and depends only on 
the ontflow loading factor. This is a very important point, 
since the relation between these two quantities is therefore a 
clear, non-degenerate tracer of the net effect of outflows (un¬ 
like the effective yield). The red lines in Fig.j^show Mz/M* 
as a function of fgas, for different values of the outflow load¬ 
ing factor. Coloured symbols in Fig. show the Mz/M* 
and fgaa observed in different radial annuli in NGC 628, 
where the colour coding represents the radial distance from 
the galaxy centre. Filled symbols are for the M08 calibra¬ 
tion, while hollow symbols are for the PP04 calibration. We 
note that this figure is just an alternative representation of 
the metal deficit discussed in the previous section (Fig. [^. 
Here we have replaced the radial distance with fgas, but the 
two quantities are related as discussed in Sec. |3.6| since the 
gas fraction increases monotonically with radius. The metal 
deficit can be visualised in Fig. as the distance between 
the observational data points and the A = 0 (solid red) line, 
which represents the Mz = y M* case. 

By using the M08 calibration the observed points are 
well reproduced with an outflow loading factor A « 1, with 
a scatter in the range 0.5 < A < 2, except for the ‘anoma¬ 
lous’ point at 3 kpc, which requires a lower loading fac¬ 
tor approaching zero. The PP04 calibration requires slightly 
higher loading factors (except again for the anomalous point 
at 3 kpc). 

We emphasise that, since we use the total mass of metals 
relative to the mass of stars, our analysis provides an esti¬ 
mate of the average outflow loading factor during the galaxy 
lifetime, and not the current outflow rate in the galaxy. In¬ 
deed, within the picture of hierarchical formation of struc¬ 
ture, it is possible that this average loading factor is mainly 
driven by the high efficiency of metal loss at higher red- 
shift, when the physical condition in the galaxy and its sub¬ 
components would have been significantly different. 


5.3 Modelling the gas metallicity 

As shown in Appendix]^ the chemical evolution framework 
adopted in the previous section also leads to a direct re¬ 
lation between gas metallicity (Zg) and gas fraction (fgas) 
that depends only on the outflow loading factor and not on 
the inflow rate and star formation efficiency. The relation 
between Zg and fgaa has been studied by several authors, in¬ 


cluding recently Ascasibar et al. (20141 and Kudritzki et al. 


(20151, who extend their analysis to resolved scales making 


use of data from NGC 628. 

However, unlike the metals locked in stars, the gas phase 
metallicity is more subject to recent gas flow events (possi¬ 
ble stochastic variations away from the ideal gas regulatory 
model), hence making it a less reliable tool for inferring the 
time averaged properties of the system. 

Keeping these caveats in mind, the model relations be¬ 
tween Zgas and fgaa are shown in Fig.[^ with the same out¬ 
flow loading factors as in Fig. The data points represent 
the observed values in NGC 628 and are colour-coded by ra¬ 
dial distance. The two different metallicity calibrations are 
identifi ed, as in Fig. by so lid (Maiolino et al. 20081 and 
hollow ( Pettini fc Pagel[2004 ( circles. The comparison of the 
models with the data points suggests outflow loading factors 


consistent with zero, or generally lower than those derived 
from the Mz/M* vs fgaa relation (Fig. |^. 

We note that the mismatch applies (with about the 
same magnitude) both to the central regions, where the to¬ 
tal mass of metals is dominated by the stellar component, 
and in the outer regions, where the total mass of metals 
is dominated by the gaseous component, therefore exclud¬ 
ing that the mismatch is simply associated with a different 
metallicity scales or calibration issues associated with the 
calculation of one of the two quantities. 

As discussed above, a possible solution can be attained 
assuming that the gas metallicity might not be representa¬ 
tive of the average evolutionary processes during the life 
of the galaxy since it can be subject to recent galactic 
‘weather’. In particular, a process that may strongly affect 
the gas metallicity, without significantly affecting the total 
content of metals, could be recent inflow of enriched gas from 
the halo. We discuss this option further in S ec. 

Finally we note that Kudritzki et al. ( 2015| have in¬ 
ferred a similarly low loading factor (A = 0.2-0.3) for NGC 
628, by fitting the gas phase metallicity gradient and the gas 
fraction. Since they use Te-based metallicities, their results 
should be compared with our analysis using the PP04 cali- 
bration,with which they are in qualitative agreement. How¬ 
ever, Kudritzki et al. (20151 make use of a chemical evolu¬ 


tion model where the inflow rate is proportional to the SFR. 
As discussed in Appendix this assumption introduces a 
dependance of the metallicity-gas fraction relation on the 
inflow rate, which is absent in our model. 


6 DISCUSSION 

In the previous section we have demonstrated that a) overall 
~ 50% of the oxygen produced by the NGC 628 out to R = 
7.0 kpc is unaccounted for and b) average outflow loading 
factors of order unity can explain the observed metal deficit, 
although there is a tension with the observed gas metallicity. 
While we cannot entirely exclude that possibility that the 
tension is due to the unavoidable systematic uncertainties in 
the determination of gas phase and stellar metallicity and 
the choice of nucleosynthetic yield/IMF, in this section we 
discuss possible ways to interpret the discrepancy in light of 
different physical mechanisms. 


6.1 Gas flows and disc-halo interaction 


It is important to remember that the chemical evolution 
framework used in this work is an intentionally simple pre¬ 
scription. In particular, while the model allows metal rich 
outflows and pristine gas inflows, and hence some exchange 
of gas and metals between the different radial bins, it does 
not treat the case of incoming metal-rich gas or changes in 
the stellar mass component due to stellar migration. 

A number of physical processes, for example, might be 
responsible for funnelling metal-rich gas into the central 
galactic regions, including viscous flows generated by gravi¬ 


tational instability or cloud - cloud collisions (Lacey & Fall 


1985 Thon fc Meusinger|1998 Ferguson fc Clarke||2001 ( or 

the presence of a bar or other non-axisymmetric perturba¬ 


tion (Minchev & Famaey 2010 Matteo et al. 20131. The 
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Figure 7. Gaseous metallicity 12 + log(0/H) versus gas fraction. 
Symbols show the values observed in radial annuli in NGC 628, 
where the colour-coding gives the galactocentric radial distance, 
as indicated in the colour bar. Solid symbols correspond to the 
mass of metals in which the ISM metallicity is calculated by using 
the M08 calibration, while open symbols use the PP04 calibration. 
The red lines correspond to the prediction of simple chemical 
evolution models with different values for the average outflow 
loading factor (A). Note that the model relations are independent 
of the inflow rate $ and of the star formation efficiency e. 


breaks seen in the stellar metallicity, and possibly the cen¬ 
tral flattening of the gas phase metallicity profile when using 
the PP04 calibration, combined with the presence of a star 
formation ring just inside the break radius, may indeed sup¬ 
port the hypothesis of an oval perturbation playing a role in 
the evolution of the disc in NGC 628. 

More importantly, enriched gas expelled by previous 
generations of supernovae may fall back onto the disc, gen¬ 
erating a galactic fountain. In the literature, however, there 
is no agreement on the details of the fountain mechanism 
(Melioli et al. 2008 2009 Spitoni et al. 20101, the radial 
distance covered by a galactic fountain, the time needed for 
the gas to rain back onto the disc and the overall impact of 
this mechanism on the chemical abundances of the disc it¬ 
self. In hydrodynamical simulations, if a hot halo is present, 
metal-enriched gas can fall back towards the galactic centre 
I Melioli et al.|2009| ). 

The presence of large amounts of metals in the halo 
of disc (and elliptical) galaxies has been inferred by recent 
HST-COS observations of absorption systems along the line 
of sight of background quasars (e.g. Tumlinson et al.]|2011 


Werk et al.|2014 1. However, these observations have probed 


the halo of galaxies at large radii (~100 kpc). In the inner 
region of the halo the gas metallicity is not known observa- 
tionally, but simulations expect that the metallicity of the 


diffuse gas should reach values higher than solar | Marinacci 
et al.|2dl4 1. 


We have attempted to model an episode of metal- 
enriched accretion at late times by modifying the analytical 
description of the metallicity evolution discussed above with 
the introduction of an accreted enriched gas mass, propor¬ 
tional to the mass of baryons already present in the galaxy, 
and with a given metallicity Zgas-accr- The blue lines in 
Fig-d shows the effect of such a late enriched accretion 
by assuming that the accreted gas is 1, 5 and 10% of the 
mass in baryons already present in the radial bin and that 
the metallicity of such enriched inflow is Zgas-accr = 3 Zq 
( which is consistent with the enrichment of the hot gas phase 
expected by some models, [Marinacci et al.]|2014[ ). Clearly 
such a model can now broadly reproduce the trend of the 
gas metallicities in NGC 628 at several radii. At the same 
time, Fig. d illustrates that such a late enriched accretion 
has little effect on the total content of metals, hence having 
little effect on the global, average outflow rate that we have 
inferred from the metal dehcit. 


In these models we have assumed an outflow loading fac¬ 
tor A = 2.0, i.e. higher than inferred in the previous section, 
to better illustrate that late enriched accretion can have a 
strong effect on the final gas metallicity but little effect on 
the total content of metals. We note that by starting with 
a lower loading factor (e.g. A = 1.0) the enriched inflow can 
reproduce the data with a smaller amount of accreted gas 
(~3-5%). 


Both gas flows in the disc and galactic fountains might 
therefore help resolve the tension between the metal deficit 
and the outflow loading factor implied from the gas phase 
metallicity for the central regions of NGC 628 by increas¬ 
ing the metallicity of the central regions beyond what is 
expected in the absence of these metal rich inflows. We note 
that there is a fundamental degeneracy in this type of model 
between the outflow loading factor and the mass (and metal¬ 
licity) of the enriched inflow when using the gas phase metal¬ 
licity alone. However better constrains can be derived if one 
models, as we do here, both the stellar and the gas phase 
metallicity at the same time, because the gas phase metal¬ 
licity is largely independent of gas flows at late times. 


areuar migration ||rtosKar er ai.|zuuoa[ 1 1 oacnim 




& Debattista 2012 Spitoni et al. 20151 might also play a role 


in alleviating the observed tension. If stars diffuse out of the 
radial bin in which they are born, we would be making incor¬ 
rect inferences for the gas fraction and oxygen mass in the 
bin. We can model this effect naively by considering that the 
observed stellar mass in each bin should be corrected for the 
net effect of migration. We experimented with this prescrip¬ 
tion and found that changes in the stellar mass up to 0.3 dex 
only generate small changes in Mz/M* for all bins within R 
= 5 kpc, since the metal budget is dominated by the stellar 
component and therefore Mz/M* is only weakly dependent 
on stellar mass. At larger radial distances the oxygen mass 
budget is dominated by the ISM and the effect of chang¬ 
ing the stellar mass content is larger. If a large fraction of 
the stellar mass currently observed at large radii is the re¬ 
sult of stellar migration from the centre, we would require a 
lower outflow loading factor in the outer disc, going towards 
a reduction of the tension between the metal dehcit and the 
observed gas metallicity. 
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Figure 8. Same as Fig. and but including a late time inflow of gas rich material with metallicity 3 Zq. The red line is the same 
as in Fig. 1^ and and represents the model with A = 2, while the blue lines represent an increasing mass of infalling metal-rich gas. 
Going upwards in both figures the blue lines represent a infalling metal-rich gas mass that is a fraction of 1%, 5% and 10% of the mass 
of baryons already present in the disc at a given location. 


6.2 Towards larger galactocentric distances 

The behaviour of the metal budget at even larger radii than 
the ones probed by our investigation constitutes a particu¬ 
larly interesting avenue for future work. Since Hi discs ex¬ 
tend substantially more than the stellar component, low lev¬ 
els of chemical enrichment in the outer disc might mean that 
a large fraction of metals have escaped our budget. Unfortu¬ 
nately there is no consensus in the literature on the nature 
of the metallicity gradient for R > 2 Rc. While |Bresolin,| 
[Kennicutt fc Ryan- Weber] ( [2012] ) and [Sanchez et ^ (20141 
argue for a flattening of the gradient beyond R = 2 Re, 


Moran et al. 120121 observed a sharp drop in the metallicity 


gradient for a subset of Hl-rich galaxies. 

Direct measurements of the chemical abundances at 
large galactocentric radii have been recently obtained for 


NGC 628 by Berg et al. (20151. Their observations contain 


12 Hll regions at radii larger than those sampled by the 
IFS data presented in this work. With the very scarce and 
scattered available data, there is no evidence for either a 
flattening or a break in the metallicity gradient of NGC 628 
out to 10 kpc. This observations are interesting in the con¬ 
text of the metal budget because they confirm the trend 
that regions at large galactocentric distances have retained 
a larger fraction of their metals compared to the inner re¬ 
gions of galaxies. In particular, relying on the metallicity 
gradient reported in [Berg et al.j ( |2015| ) we conclude that 
the radial bin around 10 kpc has retained ~80-100% of the 
metals produced in situ. 

Intriguingly, previous work has also presented evidence 


that Hll regions at extreme galactocentric distances (R > 
R24 « 3 Rc) are more metal rich than expected from in situ 
star formation (Bresolin et al.|[2009 Bresolin, Kennicutt & 


Ryan-Weber||2012 1 . These observations might be consistent 

with an extrapolation of the trend discussed above and in 
general with the idea that gas rich, relatively un-evolved 
systems have a smaller metal deficit, or might even have ac¬ 
quired metals through enriched inflows from the central re¬ 
gions (through a galactic fountain). In the context of inside- 
out disc growth, these outer regions are still in the process 
of assembly and thus represent an ideal laboratory to study 
disc formation as it happens. 


6.3 The effect of the IMF and the nucleosynthetic 
yield 


In this work we have used the [Kroupa, Tout fc Gilmor^ 
(19931 IMF, since this is the IMF which fits best the chem¬ 
ical evolution of galaxy discs ( Romano et al.|2005 20101. A 
different IMF would have an impact on our calculations by 
implying a different population-averaged net yield, since dif¬ 
ferent elements are produced by stars of different masses. In 
this section we briefly discuss the implications of adopting 
different IMFs. 

A Salpeter IMF would give a significantly higher metal 
deficit, implying a loss of about 75% of the metals produced. 
The inferred average loading factor would be about A ~ 2—4. 


In the case of a Chabrier (20031 IMF the metal deficit is even 


more severe: more than 85% of the metals should be lost and 
in the implied average loading factor should be in the range 
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A > 5. The Kroupa (20011 IMF gives results very similar to 


Chabrier 120031. We note that in going from a Kroupa to a 


Salpeter IMF we correct the stellar mass profile by a factor 
of 0.26 dex. 

Such large loading factors are extreme, especially con¬ 
sidering that these are loading factors averaged during the 
lifetime of the galaxy, and have never been seen in galac¬ 
tic winds studies locally or even in high-z galaxy discs (e.g 
Steidel et al.|[2M0 Genzel et al.|[20TT |, which instead infer 
outflow loading factors more in line with the results obtained 


by using the Kroupa, Tout & Gilmore (19931 IMF, further 


confirming that this IMF is the most appropriate for mod¬ 
elling galaxy discs. 

However, it has been pointed out that the chemical 
properties of spheroids (including bulges) may be more prop¬ 
erly described by a Chabrier ( |2003[ ) IMF. If we apply the 
latter IMF only to the central region, we get a large metal 
deficit (~ 90%) for the bulge. This would be in line with 
a similar estimates for the Milky Way, according to which 
the bulge of our Galaxy has lost about 80% of its metals 


(Greggio & Renzini 20111. The Chabrier (20031 applied to 


the central region would also imply a central outflow load¬ 
ing factor of about A ~ 5. Since we are only measuring the 
time-averaged outflow loading factor, it is not implausible to 
find that the central regions of galaxies might have experi¬ 
enced stronger outflows, despite the fact that at present time 
they sit at the bottom of the galaxy’s potential well. Indeed 
outflow rates with such large loading factor are obseerved 


in the central region of local and high-z galaxies (Feruglio 
et al.||2010[ |Rupke fc Veilleux|[2011[ [Maiolino et al.||2012[ 


Cicone et al. 


2014 


20151, and may be explained with the 


additional boost to the outflow rate associated with quasar- 
driven winds or with extreme nuclear starbursts in the past, 
when these feedback processes were likely more violent than 
in the present Universe. 

6.4 The oxygen deficit in context 

Previous work has estimated the oxygen (or metal) budget 
in star forming disc galaxies by statistical arguments, mak¬ 
ing use of scaling relations (between mass and SFR, between 
gas fraction and stellar mass, and between mass and metal- 
licity) demonstrated to hold in large samples of galaxies ob¬ 
served spectroscopically (Bouche et al. 2007 Zahid et al. 
2012[ [Peeples et al. 20141. These works predict an oxygen 
deficit ranging from 35% to 90%, with large statistical un¬ 
certainties, related both to the intrinsic scatter in the scaling 
relations used as input for the models and to the difficulty 
to robustly correct for aperture effects. Nonetheless it is re¬ 
warding to note that our estimate for the metal deficit in 
NGG 628 falls inside the range of estimates from previous 
works. 

Due to the difficulty of studying the stellar metallic- 
ity and thus deriving the total metal content in galaxies, 
the relation between gas phase metallicity and gas fraction 
has been widely used as an observationally more accessi¬ 
ble tool to investigate the impact of outflows. Recent work 


(Lilly et al. 2013 Dayal, Ferrara & Dunlop 2013 Ascasi- 
bar et al.||2014| Lu, Blanc fc Benson||2015 1 generally agrees 
on loading factors of order unity for galaxies of stellar mass 
log{Mi,/M q) ~ 10 . 

Extending this study to a larger sample of galaxies with 


resolved gas masses and metallicity information would pro¬ 
vide the logical way forward to studying the effect of out¬ 
flows across the whole galaxy population and uncover trends 
with galaxy parameters (stellar mass, SFR etc). In partic¬ 
ular, low mass galaxies (log(M*/M0) < 9) are predicted 
to exhibit higher outflow loading factors in order to match 
their low observed baryon fractions ( McGaugh et al.|2010 |, 
their metallicity ( Lu, Blanc fc Benson||20T5 1 and to recon¬ 
cile the faint-end slopes of the dark matter and galaxy mass 
functions ( Lilly et al.|2013 1. Some observational evidence for 
higher outflow rates in dwarf starbursts can be found in the 


literature (e.g Martin, Kobulnicky & Heckman 20021, but 
more evidence in this mass regime is urgently needed. 

While such a program is observationally costly, the 
synergy between large IPS surveys (CALIFA, SAMI and 
MaNGA) and radio and sub-mm interferometers (ALMA, 
NOEMA, JVLA) will greatly increase the availability of 
matched-resolution datasets to use for this type of study 
in the near future. 


6.5 A link between outflows and baryon fraction 

The metal deficit is tightly related to the well-known ‘miss¬ 
ing baryon problem’ (e.g. [Fukugita, Hogan fc Peebles|1998| 


McGaugh et al. 2010 Shull, Smith & Danforth 20121. 


Within our ‘baryon cycling’ framework we cannot account 
for the inventory of cosmic baryons, a large fraction of which 
are believed to exist in the IGM and circumgalactic medium. 
However, within the gas regulatory model, outflows modu¬ 
late the final baryonic mass of the system, while the gas 
inflow rate can be assumed to be directly proportional to 
the dark matter halo mass growth rate. [Peng fc Maioliiiol 
( |2014[ ) provide a useful parametrisation to relate the halo 
mass (Mh) and the gas inflow rate (4>) 


“h — /& fgal 


dMh 
dt ’ 


(13) 


where fb = 0.155 (Planck Collaboration et al. ][20l4| is 
the cosmic baryon fraction and fgai is the fraction of accret¬ 
ing baryons that penetrates into the inner galaxy disc, which 
could be a complex function of halo mass and redshift. Rely¬ 
ing on results from the hydrodynamical simulations in |Dekel| 
et al. (20131, we follow Peng & Maiolino (2014l and assume 


for simplicity fgai ~ 0.5 in this section. 

The evolution of the total baryonic content of the sys¬ 
tem (Mb = Mg -I- M*) is the dictated by Eq. [T||^ so that 
dMb/dt = T — 4?. Combining this with Eq. 13 we obtain 


Mb = fb fgai Mh — 


A 


1-R 


M^ 


(14) 


Rearranging we finally obtain a relation between gas 
fraction (fgas), outflow loading factor (A) and the detected 


baryon fraction {fd = 


as defined in 


20101 


McGaugh et al. 


fci = 


Mb 
fb Mh 


— fgai 1 + 


_A_ 

(1 — R)(l -I- fgas) 


(15) 


Using representative values for NGG 628, fgas = 0.2 
and A = 1 we obtain fd ~ 0.23. This very rough order-of- 
magnitude calculation leads to a detected baryon fraction 
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16 F. Belfiore et al. 


that is surprisingly close to the detected baryons fractions 
measured via the Tully-Fisher relation in |McGaugh et ah] 
120101, who quote/d ~ 0.16 —0.19 for disc galaxies of stellar 
mass log(M*/M0) = 9.9 — 10.1, comparable to NGC 628. 


7 CONCLUSIONS 


In this work we have presented a detailed study of the metal 
content and distribution in the nearby disc galaxy NGC 628. 
We have mapped the gas phase metallicity by exploiting the 
largest integral field spectroscopy mosaic available to date, 
thus enabling us to extend our study out to a galactocentric 
radius of 7 kpc (~ 3 Rc). We have combined the gas metal¬ 
licity map with the stellar metallicity radial prohle presented 
in previous work ( Sanchez-Blazquez et al.|2014 |, the stellar 
mass surface density distribution E* (inferred from exten¬ 
sive multi-band photometry) and the observed gas surface 
density (Ehi,Eh 2) inferred from CO and Hi maps. 

By comparing the mass of metals observed in the 
gaseous and stellar components with the mass of metals pro¬ 
duced by stars in the same region (inferred assuming the 
latest determination of the oxygen nuncleosynthetic yield, 
Vincenzo et al.|2015 1 we have obtained a detailed, spatially 
resolved metal budget as function of galactocentric radius. 
The main results from this analysis are the following: 

(i) On average about 50% of the metals produced by the 
stars in each galactic region have been lost. The fraction of 
metals lost is higher (~70%) in the central region, which is 
dominated by the bulge. There is also tentative indication 
that the fraction of metals lost decreases at large galactocen¬ 
tric radii. Cumulatively, out to a radius of 7 kpc, 45%-50% 
of the total amount of metals produced have been lost by the 
galaxy. Evidence from recent measurements of outlying Hll 
regions from Berg et al. ( 2015| ) (out to 10 kpc) confirms the 
trend that regions in the outer disc have retained a higher 
fraction of the metals they have produced, with regions at 
around 10 kpc retaining ~ 80-100 % of the metals produced 
there. 

(ii) We have used simple (‘bathtub’ and ‘gas regulatory’) 
analytical models involving gas outflow and inflow to model 
the data. In this framework the relation between metals- 
to-stellar mass ratio (Mz/Mgtar) and gas fraction (fgaa) is 
independent of the inflow rate and of the efficiency of star 
formation, and depends only on the outflow loading factor 
in models with constant inflow rate. We find that an out- 
fiow loading factor A « 1 can explain the data at most radii, 
though with large scatter (0.5 < A < 2). We emphasise 
that, since this inference is based on the total metal bud¬ 
get, we derive a time-averaged outflow rate, representative 
of the strength of outflows over the lifetime of the galaxy. 
This information is thus complementary to that obtained 
in studies based on the observed gas kinematics, which are 
sensitive only to the ongoing/recent outflow activity. 

(iii) A larger loading factor may be derived if we make 
use of the Chabrier IMF. This IMF might be relevant for 
the central spheroid, thus implying a higher outflow rate 
during the early formation of the bulge. 

(iv) The gas phase metallicity is more sensitive to recent 
evolutionary processes because it is associated with a com¬ 
ponent of the galaxy (the ISM) that is potentially subject 


to rapid evolution. We hnd that the observed gas metallic¬ 
ity cannot be reproduced by the same parameters inferred 
from the metal budget, independently from the choice of 
metallicity calibrator used. In particular, the observed gas 
phase metallicities would require little or no outflow. We 
show that the apparent tension could be explained with in¬ 
fall of enriched gas from the halo at late times. Indeed, such 
a process would strongly affect the observed gas metallicity, 
but not the total mass of metals. 

(v) Finally, we warn against the use of analytical models 
in which the inflow rate is assumed to be proportional to 
the star formation rate. This assumption, adopted in pre¬ 
vious work, is not physically motivated and introduces an 
unnecessary degeneracy between inflow and outflow rates 
when modelling chemical abundances. Analytical solutions 
are presented in this work for models assuming a constant, 
or slowly varying, accretion rate, which is a much better 
representation of gas inflows in the cosmological framework. 
However, further work is needed to reliably test the valid¬ 
ity of the underlying assumptions of the simple regulatory 
model presented in this work for describing regions within 
galaxies. 
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APPENDIX A: ANALYTICAL SOLUTION FOR 
CONSTANT INFLOW MODELS 


In this section we summarise the chemical evolution formal¬ 
ism of |Peng fc Maioliii^ ( |2014[ ) and present further exact 
solutions for the metallicity and total metal content. 


Peng & Maiolino (20141 models offer solutions for equa¬ 
tions Eq. under the following assumptions: 


(i) Constant inflow rate $. 

(ii) A linear relation star formation law SFR = e Mgas, 
with constant e. 

(iii) Outflow rate proportional to the SFR through a con¬ 
stant loading factor A: 4^ = A SFR. 


We further assume that the net nucleosynthetic yield p 
and return factor R are constant with time and the the in¬ 
flow is of primordial gas (with zero metallicity). Under these 
assumptions analytical solutions can be obtained with an 
explicit time dependance. As discussed in |Peng fc Maiolino| 
(20141 the equilibrium timescale 


^“£(l-R-fA) 

is the natural timescale driving the chemical evolution 
of the system. For t >> r physical properties of galaxies 
tend to their ‘equilibrium’ values. As discussed in |Peng fc| 


Maiolino (20141, low mass dwarf galaxies and chemically un¬ 


evolved systems will not satisfy the equilibrium condition 
and need to be studied in the t < r regime. 

In Table [XT] we summarise the analytical solutions for a 
number of galaxy properties. We note that, apart from the 
assumptions stated above, the solutions are mathematically 
exacl[^ We also remark that, as expected, the only timescale 
responsible for the chemical evolution of the system is the 
equilibrium timescale dehned above (Eq. 

The limits of the solutions for the equilibrium case 
(t >> r) have a simple physical interpretation. Star for¬ 
mation exactly balances out the inflow and outflow rates, 
the gas mass (and hence the SFR) stays constant and the 
stellar mass grows linearly with time. In this regime the gas 
phase metallicity plateaus at a constant value 


® The solution to the gas phase metallicity (Zg) presented here 
differs from the one in |Peng &: Maiolino| ( |20T4| |, who solve for the 
metallicity evolution making further approximations to simplify 
the algebra. 


Zg{eqm) ^per = ^ (A2) 

which depends only on the net yield, the return fraction 
and the outflow loading factor. Note also that the equilib¬ 
rium value of the Mz/M* ratio is (1 — £ A r), lower than 
the expected value Mz/M* = in absence of outflows 
(i.e. for A = 0). The gas fraction fgas can be considered as 
a function of t,r and £, and is the best property to quan¬ 
tify the degree of chemical evolution of the system: systems 
with low gas fraction have not reached equilibrium and are 
chemically un-evolved, while systems with low gas fraction 
are in equilibrium and chemically evolved. 

The analytical solutions are also an important tool to 
study the parameter dependences of different physical prop¬ 
erties. In particular we can easily show that the relation 
between Mz/M* and fgas depends only on the value of the 
outflow loading factor, and does not depend on the assumed 
inflow rate and star formation efficiency. In other words the 
position of a point in the fgas vs Mz/M* plane depends only 
on A. To prove this it is sufficient to observe that both fgas 
and Mz/M* are only a function of A and t/r. Therefore, 
eliminating t/r, it would be possible to write a relation be¬ 
tween fgas and Mz/M* which depends only on A. Unfortu¬ 
nately the complex nature of the solutions does not allow 
such relation to be written in closed form. Physically this 
statement can be interpreted as stating that the metal deficit 
(i.e. Mz/M*) depends only on how evolved the system is (i.e. 
its fgas) and the strength of its outflows (i.e. A). As shown in 
App. 0 this simple result is not valid anymore if we assume 
that the inflow rate is proportional to the SFR. 

Finally we note that exactly the same argument can 
be applied to the relation between fgas and Zg, which can 
therefore be shown to depend only on the outflow loading 
factor and not on the inflow rate and the star formation 
efficiency. 


APPENDIX B: EQUILIBRIUM, TIMESCALES 
AND INFLOW RATES IN DIFFERENT 
‘BATHTUB’ MODELS 


In this section we briefly recap the main features and differ¬ 
ences between the growing set of chemical evolution ‘bath- 


tub’ models. We will focus on a comparison of the models 

presented in Dave, Finlator & Oppenheimer 

(20111, Lilly 

et al.| (|2013P, |Dayal, Ferrara & Dunlop] (|2013 

and |Peng & 


Maiolino (20141, the latter being our reference model. The 


casual reader is encouraged to read the summary at the end 
of this section. 

A ‘bathtub’, or ‘gas regulatory’ model, is defined here 
to be a system assumed to obey the equations [TJ|^ An ‘ideal’ 
bathtub model is one in which we assume 


(i) A linear relation star formation law SFR = e Mgas, 
with constant e. 

(ii) Outflow rate proportional to the SFR through a con¬ 
stant loading factor A: <& = A SFR. 

the one presented in 
), is an ‘equilibrium’ 
model. An equilibrium model is dehned by the assumption 
that the gas mass in the system does not change with time 


The simplest bathtub model, like 


Dave, Finlator & Oppenheimer (2011 
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Table Al. Exact analytical time-dependent solution of t he|Peng &: Maiolino| | [2014[ | models with constant inflow rate for different galaxy 
properties. The equilibrium timescale r is defined in Eq. 1^3 


Galaxy property 


Exact solution 


t << r 


t >> r 


M, 

SFR = e Mg 
M* 

f = ^3 

Jgas — Mg+M* 

Zg = Mzg/Mg 

Mz 

Mz/M^, 


<!> T (1 - 

£ <I> r (1 - 

$ £ (1 - R) (t/r - (1 - e-‘/^)) 

1 + £t(1-R) U — 


p £ T 1 - 


t/r e 


-t/7 


\-&-t/r 


P 

1-R 


1 — £ A T I 1 


((l+t/r) e 

t/r+e-t/r-l 


$ t 
£ $ i 

<*>£ (1-R) *2/2 


7- 

£ $ T 

<I> £ (1 — R) T i 


(l + £ (1-R) i/2)-i (l + £ (1-R) t)-i 
p £ t/2 


> £ <!> t2 [(1 - e A t)(t/T - (1 - e-‘/^)) - £ A t{{t/T + l)e-*/^ - 1)] 


p £ <I> *2/2 


P 

1-R 


per 

pe^{l — £\t)Tt 


(dMgas/dt = 0). This assumption is justified by Dave, Finla- 


tor & Oppenheimer (20111 on grounds of its relevance to cos¬ 


mological hydrodynamical simulations. In equilibrium Eq. 
reduces to 


$ = (1 - i? -f A) SFR (in eqm). 


(Bl) 


Note that even in an equilibrium model the SFR will 
change in time, however equilibrium requires that it does so 
at the same rate as the inflow rate. The same conclusion is 
evident in the equilibrium solution for the SFR presented in 
Table lAl] above. 


Dave, Finlator & Oppenheimer (20111 also argue that, 


when departures from equilibrium are allowed, the bathtub 
model naturally generates a secondary dependance of the 
mass-metallicity relation on SFR, in the sense that more gas 
rich systems will have higher SFR and lower metallicities. 
This feature of the model is generally considered a success 


(Lilly et al. 20131, since it is capable of reproducing the 
observed fundamental metallicity relation ([Mannucci et al.] 


2010 Andrews & Martini 20131. Moreover this feature is 


mostly independent of the detailed time-evolution of the cos¬ 
mological inflow rate, since it is just a general characteristic 
of the non-equilibrinm bathtub models. 


Subsequent work in Lilly et al. (20131 has shown that 


equilibrium is generally a good assumption for gas-poor, 
chemically evolved, massive low-redshift galaxies, while it 
is not a good assumption for the gas-rich dwarfs or galaxies 


at high redshift (see in particular Sec. 2.3 and 2.4 of Lilly 
let al.l20T3l for a detailed discussion of the relevant timescales 
and Sec. 5 of [Peng fc Maiolino||2014| for an estimate of the 
equilibrium timescale for gas rich systems). It is therefore 
clear that in applying the bathtub model to the chemically 
un-evolved, gas rich outer disc of NGC 628 a model that 
takes into account departure from equilibrium is needed. 

The simplicity of the framework offered by the bathtub 
model has encouraged several authors to seek analytical so¬ 
lutions for its fundamental equations. Unfortunately, if we 
do not assume the system is in equilibrium, the relation be¬ 
tween metallicity and gas fraction (which can be interpreted 


as natural ‘clock’ in the context of chemical evolution) will 
depend on the assnmed time evolution of the inflow rate. 

A notable special case is represented by the popular 
solutions for the case in which in the inflow rate is propor¬ 
tional to the SFR (Matteucci fc Chiosi]1983' Dayal, Ferrara 


fc Dunlop[2013 Kudritzki et al.|2015 1. The success of these 

solutions is to be attributed to the fact that they are ar¬ 
guably the simplest example of non-equilibrium models (i.e. 
dMgas/dt 7^ 0) which have closed-form analytical solutions. 
As all non-equilibrium bathtub models, they show the usual 
behaviour of tending to equilibrium at late times and pre¬ 
dict a secondary relation on SFR in the mass-metallicity 
relation ( Dayal, Ferrara fc Dunlop|2013 1. However, we warn 
the reader against fallaciously reversing the argument, i.e. 
the existence of the fundamental metallicity relation does 
not necessarily imply that the inflow rate must be propor¬ 
tional to the SFR! 

It is important to stress that, while the other assump¬ 
tions made in the context of the bathtub model can be 
justified, albeit approximately, on physical or observational 
grounds, there is no physical or observational reason to be¬ 
lieve that the inflow rate should follow the same time evo¬ 
lution as the SFR when the system is not in equilibrium. 
Confusion often arises in the literature because, in equilib¬ 
rium, all bathtub models predict the time evolution of the 
SFR to the follow that of the inflow rate. However, this con¬ 
clusion is only true when the system is at equilibrium and 
not trne in general as the system approaches eqnilibrium. 

On the other hand, the cosmological inflow rate follows 
a well-defined trend with halo mass and redshift. Both semi¬ 
analytics and cosmological simulations demonstrate that the 
halo accretion rate can be parametrised as 


dt 




(B2) 


where Mh is the halo mass and a ~ 0.030 Gyr and 


P ~ 2.5 (Dekel et al. 20131. For an Einstein-de-Sitter uni¬ 
verse the Eq. |B2| admits a particularly simple solution lead¬ 
ing to 
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Mf = a = a/{Hon]fl) (B3) 


where the subscripts f and and i refer to final/initial 
ma sses and redshifts. Therefore, by substituting back into 
Eq. B2 the halo mass accretion rate at redshift z, Mh{z) is 


given by 


Mh{z) = 0.03 Mi + zf 


(B4) 


As reported in previous work (Dekel et al.||2013 Lilly 


et al. 20131, this halo accretion rate does not evolve very 


much over most of the age of the Universe. In particular, 
it changes by about a factor of 2 between z = 5 and 0.3. 
We neglect here that fact than only a fraction of accret¬ 
ing baryons will reach the galaxy disc and are thus avail¬ 
able for star formation, and this fraction might be a com¬ 
plex function of mass, redshift and perhaps other quantities 
(environment). Under these assumptions, |Peng fc Maiolino| 
(20141 argue that the inflow timescale is much longer than 


all other timescales relevant to galaxy evolution (including 
the equilibrium timescale) and, therefore, it is an excellent 
approximation to assume the gas inflow rate to be constant 
in time. By making the assumption that the timescale for 
changes in the infall rate is much longer than the equilib¬ 
rium timescale we obtain the models presented in Table [Alj 
which have been used throughout this paper. 

To demonstrate the physical difference between the so¬ 
lutions that assume a constant inflow rate and the solutions 
which assume the inflow rate to the proportional to the SFR 
at all times, we present a comparison of their time evolution. 
The first thing to appreciate is that, if we assume the inflow 
rate to the proportional to the SFR ($ = tu SFR), Eq. 
only admits the following solution for the gas mass Mg 


= {u-X-l + R)eMg ^Mg = 

^ (B5) 

where Mq is the initial gas reservoir at t=0. Physical 
solutions require 1 — R — u)-|-A > 0, and hence imply an 
exponentially declining SFR. The gas equilibrium timescale 
(r') for this model depends on uj and is given by 


1 

£(1 — i? -I- A — a;) ’ 


(B6) 


This timescale tends to the timescale in lAll for small w. 
We therefore refer to the model with SFR proportional to 
the inflow rate as an cu-model. The equilibrium behaviour 
of the cj-model is very different from the constant inflow 
regulator. In the oj-model equilibrium can only be reached 
for SFR = e Mg = 0. Because of its defining assumption, 
this model implicitly assumes the inflow rate ($ = cu SFR) 
to be exponentially declining (unless we wish to build a 
more general model with a different star formation law, i.e. 
SFR ^ e Mg, which we will not discuss here). In light of 
the discussion on the cosmological accretion rate above, this 
is a much worse parametrisation than the assumption of a 
constant inflow rate. However, an exponentially declining in¬ 
flow rate is still a useful model in the context of the overall 
evolution of galaxies, since it can represent the shutdown of 
star formation due to a declining inflow rate. This model is 


therefore similar to the ‘strangulation’ models used in|Peng,| 


Maiolino & Cochrane (20151 to study the build up of the 


red sequence. 

Although the relation between metallicity and gas mass 
in the cu-model has been widely used in the literature, it is 
instructive to present here the solutions for the time evolu¬ 
tion of this model, which can be found in Table [Blj We note 
that, unlike in the constant inflow model, the metallicity in 
the w-model does not come to equilibrium on the gas equi¬ 
librium timescale defined in IBGI but on the dilution timescale 
Tdii = (e = Mg/^. Therefore the final relation between 
metallicity and gas content will depend on combination of 
the gas equilibrium and the dilution timescales. In particu¬ 
lar, by combining the solutions for Zg and Mg we obtain the 
well-known relation 


^9 = 


P 

UJ 



XjUa) 


(B7) 


and by combing the solution for M* and Mg we can 
eliminate Mo from the above the re-write the metallicity in 
the cu-model completely in term of observables. 




P_ 

UJ 




(1 -I- w) 
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where w = (X — cu)/(l — R). As can be explicitly seen 
from this equation, the relation between Zg and gas fraction 
depends on the inflow rate (through to), causing the observed 
degeneracy in the determination of cv and X when using this 
model. 

In summary 


(i) Equilibrium (defined by the requirement dMgas/dt = 
0) bathtub models predict the SFR to be proportional to 
the inflow rate. 

(ii) Non-equilibrium bathtub models tend to equilibrium 
at late times. These models can naturally reproduce the fun¬ 
damental metallicity relation (i.e. imply a SFR dependance 
in the mass-metallicity relation). 

(iii) The cosmological accretion rate, as calculated using 
semianalytics or cosmological simulations, is slowly varying 
over most of the age of the Universe. Hence a constant gas 
inflow rate is generally a good assumption for the purpose 
of chemical evolution modelling. In a constant inflow model, 
the relation between metallicity and gas fraction depends 
only on the outflow loading factor and not on the inflow 
rate. 

(iv) Non-equilibrium models which fix the inflow rate to 
be proportional to the SFR implicitly assume the inflow rate 
to be exponentially declining. This is not a good represen¬ 
tation of the cosmological accretion rate for ‘main sequence’ 
star forming galaxies, but could constitute a good model 
for strangulation, as in Peng, Maiolino & Cochrane (20151. 
Moreover in this model the relation between the metallicity 
and the gas fraction depends on both the outflow loading 
factor and the proportionality factor between inflow rate 
and SFR, causing a degeneracy between inflows and out¬ 
flows when interpreting metallicity data. 
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Table Bl. Exact analytical time-dependent solution of the ideal regulator with inflow rate proportional to the SFR ($ = a; SFR) for 
different galaxy properties. Note that Mq below refers to the gas mass at t=0 and r' is the gas equilibrium timescale defined in Eq. |B6| 


Galaxy property 

Exact solution 

t = 0 

t >> r' 

Mg 

Mo 

Mo 

0 

SFR = e Mg 

e Mo 

e Mo 

0 

M* 

Mo t' 6 (l-R) 

0 

Mo r' £ (1 - R) 

r — 

Jgas — Mg+K'U 

[l + er' (l-R) -1)]"^ 

OO 

0 

Zg = Mzg/Mg 


0 

£ P 
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